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ABSTRACT 

We provide a first look at the results of the Herschel Gould Belt survey toward the IC5146 molecular cloud and present a preliminary analysis 
of the filamentary structure in this region. The column density map, derived from our 70-500 fim Herschel data, reveals a complex network of 
filaments, and confirms that these filaments are the main birth sites of prestellar cores. We analyze the column density profiles of 27 filaments and 
show that the underlying radial density profiles fall off as r _1 5 to r -2 5 at large radii. Our main result is that the filaments seem to be characterized 
by a narrow distribution of widths having a median value of 0.10 + 0.03 pc, which is in stark contrast to a much broader distribution of central 
Jeans lengths. This characteristic width of ~ 0.1 pc corresponds to within a factor of ~ 2 to the sonic scale below which interstellar turbulence 
becomes subsonic in diffuse gas, supporting the argument that the filaments may form as a result of the dissipation of large-scale turbulence. 
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1. Introduction 

Understanding how stars form out of the diffuse interstellar 
medium (ISM) on both global and local scales is a fundamen- 
tal open probl em in contemporary astrophysics (see McKee & 
|Ostr iker 2007 for a recent review). Much observational progress 
is being ma de on this front than ks to the Herschel Space 
Observatory (Pilbra tt et al. 2010| ). In particular, the first re- 
sults from the Gould Belt and Hi-GAL imaging surveys with 
Herschel have revealed a profusion of parsec-scale filaments in 
Galactic molecular clouds and suggested an intimate connection 
between the filamentary structure of the ISM and the forma- 
tion proces s of dense cloud cores ( Andre et al. 2010; Molinari 
et al.|2010| ). Remarkably, filaments are omnipresent even in un- 



bound, non-star-forming complexes such as the Polaris translu 
cent c loud ([M en'shchi kov et al.|2010[ [Miville-Deschenes et al. 
|2010[ | Ward-T hompso net al. |20 1 0| )TFurthermore, in active star- 
forming regions such as the Aquila Rift cloud, most of the 
prestellar cores identified with Herschel are located within grav- 
itationally unstable filame nts for which th e mass per unit length 
exceeds the critical value ( |Q striker || 1 9 64| ) , Mii ne?cr i t = 2c 2 s /G ~ 
15M /pc, where c s ~ 0.2 km/s is the isothermal sound speed 
for T ~ 10 K. The early findings from Herschel led | Andre 
|et al. ( 2010| ) to favour a scenario according to which core for- 
mation occurs in two main steps. First, large-scale magneto- 
hydrodynamic (MHD) turbulence generates a whole network of 
filaments in the ISM (cf. |Padoan et al.|2001[|Balsara et al.|200l) ; 
second, the densest filaments fragment into prestellar cores by 
gravitational instability (cf. Inutsuka & Miyam a[l997| ). 

To refine this observationally-driven scenario of core forma- 
tion and gain insight into the origin of the filamentary struc- 
ture, an important step is to characterize the detailed physical 



* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



properties of the filaments seen in the Herschel images. Here, 
we present new results from the Gould Belt survey obtained to- 
ward the IC5146 molecular cloud and analyze the radial den- 
sity profiles of the numerous filaments identified in this cloud. 
Based on a comparison with a similar analysis for the Aquila and 
Polaris regions, we show that the filaments of IC5146, Aquila, 
and Polaris are all characterized by a typical width of ~ 0.1 pc 
and discuss possible physical implications of this finding. 

IC5146 is a st ar-forming cloud in Cygnus located at a dis- 
tance of ~ 460 pc ( |Lada et al.|1999| but see Appendix A), which 
consists of the Cocoon Nebula, an H n region illuminated by the 
B0 V star BD + 46° 3474 , and two paral lel filamentary stream- 
ers extending to the west ( [Lada et al.|1994| ). In addition to these 
two streamers, the Herschel images reveal a whole network of 
filaments, which are the focus of the present letter. 



2. Observations and data reduction 

Our observations of IC5146 were made on 29 May 2010 in the 
parallel scan-map mod e of Herschel. An area of -1.6 deg 2 was 
covered by both PACS (|Poglitsch et al.|2010| ) at 70 /mi, 160 /mi 
and SPIRE ( [Griffin et al.|2010| ) at 250 //m, 500 /mi, with 

a scanning speed of 60 s _i . The PACS data were reduced with 
HIPE version 3.0.1528. Standard steps of the default pipeline 
were applied, starting from the raw data, taking special care 
of the deglitching and high-pass filtering steps. The final maps 
were created using the photProject task. The SPIRE observations 
were reduced using HIPE version 3.0.1484. The pipeline scripts 
were modified to include data taken during the turnaround of the 
telescope. A median baseline was applied and the 'naive' map- 
making method was used. 

Thanks to their unprecedented spatial dynamic range in the 
submillimeter regime, the Herschel images (see, e.g., online 
Fig. [5] and Fig. [6]) provide a wealth of detailed quantitative in- 
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Fig. 1. (a) Column density map derived from our SPIRE/PACS observations of IC5146. The resolution corresponds to the 
36.9" HPBW resolution of SPIRE at 500 fim. The positions of the YSQs and starless cores detected with Herschel (using get- 



sources, a source-finding algorithm described in |Men'shchikov et aL 2010b are plotted as red stars and blue triangles respectively; 
bound starless cores are in dark blue (cf. Fig. [6}) and Konyves et al. 2010| for classification details). The locations of the two fil- 
aments (6 and 14) whose radial profiles are shown in Fig. [2] and online Fig. [4] are marked by the white rectangles, (b) Curvelet 
component of the column density map, with the network of 27 identified filaments shown in blue. 



formation on the large-scale, filamentary structure of the cloud. 
It is this filamentary structure that we discuss in the following. 

3. Analysis of the filamentary structure 

3.1. Identification of a network of filaments 

Based on the Herschel images at five wavelengths, a dust temper- 
ature (7d) map (see online Fig. [6^) and a column density (Nu 2 ) 
map (Fig.[lJ) were constructeafor IC5146, assuming the dust 
opacity law of Hildebrand (1983, see Appendix A) and follow- 
ing the same procedure as [Konyves et al.| ( |2010| ) for Aquila. We 
decomposed the column density map on curvelets and wavelets 
(e.g. Starck et al.|2003| ). The curvelet component image (Fig. [I])) 
provides a high-contrast view of the filaments (after subtraction 
of dense cores, which are contained in the wavelet component). 
We then applied the DisPerSE algorithm ( Sousbie 2011) to the 
curvelet image to take a census of the filaments and trace their 
ridges. DisPerSE is a general method to identify structures such 
as filaments and voids in astrophysical data sets (e.g. gridded 
maps). The method, based on principles of computational topol- 
ogy, traces filaments by connecting saddle points to maxima 
with integral lines. From the filaments identified with DisPerSE 
above a 'persistence' threshold of 5 x 10 20 cm~ 2 in column den- 
sity (~ 5cr in the curvelet map), we built a mask on the same grid 
as the input ma p, with values o f 1 along the filament ridges and 
elsewhere (cf. |Sousbie|2011 for a definition of 'persistence'). 
Based on the skeleton mask obtained for IC5146, we identified 
and numbered a total of 27 filaments (shown in blue in Fig.[T])). 

3.2. Radial density profiles of the filaments 

To construct the mean radial density profile of each filament, the 
tangential direction to the filament's ridge was first computed at 
each pixel position along the filament. Using the original column 
density map (cf. Fig.[T^), we then measured a radial column den- 
sity profile perpendicular to the ridge at each position. Finally, 
we derived the mean radial profile by averaging all of the profiles 
along the filament (see Fig. [2]). In order to characterize each ob- 
served (column) density profile, we adopted an idealized model 



of a cylindrical filament with radial density and column density 
profiles (as a function of cylindrical radius r) of the form: 



P P (r) 



[l + (r/tf flat ) 2 ] 



P/2 



Z,(r) = A P 



Pc^fla 



[l + (r/tf fla 2 ]" 



(1) 



where X = /dmuNu 2 , H = 2.33 is the mean molecular mass, 

A p = I-oo (i+^V 2 is a finite constant factor for P > 
p c is the density at the center of the filament, and is the 
characteristic radius of the flat inner portion of the density pro- 
fil^] For simplicity, the inclination angle i of the model fil- 
ament to the plane of the sky is assumed here to be / = 
(but see online Appendix [A] for the effect of i ± 0). At large 
radii (r » R^), the model density profile approaches a power 
law: p p ( r) ~ p c (r/R^ t )~ p . An important special case is the 



|Ostriker| ( |1964 ) model of an isothermal filament in hydrostatic 



equilibrium, for which p = 4, A p = n/2, and R^ t corresponds 
to the thermal Jeans length at the center of the filament, i.e., 
#flat^j(r = 0) = c 2 /(Gi: r =o). 

When fitting such a model profile to the observations, p c , 
Rfia, and p, along with the peak position of the filament, were 
treated as free parameters. The results (see, e.g., Fig. [2] and on- 
line Table 1) show that the observed profiles are generally well 
fitted with 1.5 < p < 2.5. None of t he observed fil aments has 
the steep p = 4 density profile of the Ostriker (19641 model. A 
similar conclusion was already reported by |Lada etaL (1999]) in 
the case of the main streamer (i.e., Filament 6 here). 

The above radial profile analysis can also be used to de- 
rive accurate masses per unit length, Mn ne , for the filaments 
by integrating column density over radius: Mn m = J £ b s (r)dr. 
The main source of uncertainty lies in the difficulty of defining 
the edges of the filaments, especially in crowded regions. The 
method is nevertheless quite robust when the bounds of inte- 
gration are reasonably well defined (cf. Fig. [2]). The values of 



1 Model profiles of th e form described by Eq. (1) are sometimes 
called Plummer-like (cf. Nutter et al. 2008 and references therein). 
When /?flat is much smaller than the spatial resolution of the observa- 
tions, such profiles are effectively equivalent to power laws. 
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Fig. 2. (a) Mean radial column density profile perpendicular to the supercritical filament labeled 6 in Fig.[T|) (black curve). The area 
in yellow shows the dispersion of the radial profile along the filament. The inner curve in light blue corresponds to the effective 
36.9" HPBW resolution (0.08 pc at 460 pc) of the column density map of Fig. [T^ used to construct the profile. The dotted red 
curve shows the best-fit model of the form expressed by Eq. (1), with p = 2. For comparison, the dashed blue curve shows a model 
with p = 4, corresponding to a hydrostatic isothermal equilibrium filament (Ostriker 1964 ). The two dashed vertical segments at 
~ ±1.5 pc mark the bounds within which the profile was integrated to estimate the mass per unit length Mn nQ . (b) Mean radial dust 
temperature profile measured perpendicular to filament 6. (c) Same as (a) displayed in log-log format to better visualize the flatter 
inner part of the profile and the power-law behavior at large radii. 



Mii ne derived with this method (see online Table 1) are typi- 
cally 20% higher than simpler estimates assuming that the fila- 
ments have Gauss ian column density profiles (cf. Appendix of 
And re et "aL||2010| ). Most bound prestellar cores appear to be 
located within supercritical, gravitationally unstable filaments 
with Mii ne > Mii ne , C rit (e.g. Filament 6 - see online Fig. |6fr ), a 
similar result to that already obtained by Andr e et aL| ([20TD ) in 
Aquila. Note that both methods of estimating Mii ne (profile in- 
tegration or Gaussian approximation) yield essentially the same 
conclusion on the gravitational instability of the filaments. 

4. A characteristic width for interstellar filaments ? 

To construct a reliable distribution of filament widths (cf . online 
Fig. |7^), we applied Gaussian fits to the radial column density 
profiles, as these tend to be more robust than the more sophisti- 
cated fits discussed in Sect. 3.2. A Gaussian fit to a model pro- 
file of the form of Eq. (1) with p = 2 indicates that the derived 
FWHM width is equivalent to ~ 1.5 x (27?fl at ). It can be seen 
in Fig. [7^ that the filaments of IC5146 have a narrow distribu- 
tion of deconvolved FWHM widths centered around a typical 
(median) value of 0.10 ± 0.03 pc (see also Fig. [3]). Note that the 
same filaments span more than an order of magnitude in central 
column density (cf. online Table 1), implying a distribution of 
central Jeans lengths [Aj (r = 0) = c 2 / (GE r =o)l f rom ~ 0-02 pc 
to > 0.3 pc, which is much broader than the observed distribu- 
tion of widths (see Fig.|3]and blue dashed line in Fig. [7^). 

In order to check that the measured filament widths were not 
strongly affected by the finite resolution of the column density 
map, we performed the same radial profile analysis using the 
SPIRE 250 pm map, which has a factor ~ 2 better resolution 
(18.1" HPBW). The resulting distribution of filament widths is 
shown as a dotted histogram in online Fig. [7^ and is statistically 
very similar to the original distribution at 36.9" resolution (at the 
50% confidence level according to a Kolmogorov-Smirnov test). 

We also performed a similar analysis of the filamentary 
structure in two other regions located at different distances, the 



Aquila and Polaris fields, also observed by us with Herschel. 
The Aquila regio n is a very active s tar-forming complex at 
d ~ 260 pc (e.g. Bontemps et al. 2010, and references therein), 
while the Polaris field is a high-latitude translucent cloud with 
little to no star formation at d ~ 150 pc (e.g. |WarrpTh ompson 
et al. 2010, and references therein). Following the same proce- 
dure as in Sect. 3.1, we identified 32 filaments in Aquila and 
31 filaments in Polaris spanning three orders of magnitude in 
central column density from ~ 10 20 cm -2 for the most tenuous 
filaments of Polaris to ~ 10 23 
Aquila (see also Fig. 1 of And re et aL]|2010| ) 



2 for the densest filaments of 
The distribution 



of deconvolved FWHM widths for the combined sample of 90 
filaments is shown in online Fig.|7]3. This combined distribution 
is sharply peaked at 0.10 ± 0.03 pc, which strengthens the trend 
noted above for IC5146. While further tests would be required 
to fully investigate potential biases (especially given uncertain- 
ties in cloud distances - see Appendix A) and reach a definitive 
conclusion, our present findings suggest that most interstellar fil- 
aments may share a similar characteristic width of ~ 0.1 pc. 



5. Discussion 

The results presented in Sect. 3 and Sect. 4 can be used to discuss 
the formation and evolution of the observed filamentary struc- 
ture. Three broad classes of models have been proposed in the 
literature to account for filaments in molecular clouds, depend- 
ing on whe ther global gravity (e.g . Heitsch et al.|2008| ), magnetic 
fields (e. g. Nakamura & Li 2008), or large -scale turbulence (e.g. 
|Padoan et al.||2001[|Mac Low & Klessen||2004| ) play the domi- 
nant role. Large-scale gravity can hardly be invoked to form fila- 
ments in gravita tionally unbound complexes s uch as the Polaris 
flare cloud (cf. Heithause n & Thaddeus|1990| ). While magnetic 
fields may play an important role in channeling mass accumu- 
lation onto the densest filaments (e.g. Goldsmith et al. 2008, 
Nakamura & Li 2008), our Herschel findings appear to be con- 
sistent with the turb ulent picture. In the scenario proposed by 
Padoan et al. ( 2001| ), the filaments correspond to dense, post- 
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shock stagnation gas associated with regions of converging su- 
personic flows. One merit of this scenario is that it provides an 
explanation for the typical ~ 0.1 pc width of the filaments as 
measured with Herschel (Sect. 4 and Fig. [3}. In a plane-parallel 
shock, the thickness A of the postshock gas layer/filament is re- 
lated to the thickness L of the preshock gas by A « L/M(L) 2 , 
where M(L) is the Mach number and M(L) 2 is the compression 
ratio of the shock for a roughly isothermal radiative hydrody- 
namic shock. Thus, the thickness A is independent of the scale L 
given Larson's linewidth-size relation [Ai(L) oc <x v (L) oc L 5 ]. 
In this picture, the postshock thickness of the filaments effec- 
tively corresponds to the sonic scale A s at which the 3D turbu- 
lent velocity dispersion equals the sound speed (i.e., M(A S ) = 1), 
leading to A « A s ~ 0.05-0.15 pc according to recent determina- 
tions of the linewidth-size relationship in molecular clouds (e.g. 



Heyer et al. 2009; Federrath et al.|2010|). 



If large-scale turbulence provides a plausible mechanism for 
forming the filaments, the fact that p restellar cores ten d to form 
in gravitationally unstable filaments ( An dre et al.|2 010) suggests 
that gravity is a major driver in the subsequent evolution of the 
filaments. The power-law shape of the outer density profiles, 
especially in the case of supercritical filaments (Sect. 3.2 and 
Fig. [2]), is also suggestive of the role of gravity. Although the ob- 
served profiles are shallower than t he profile of a self-gravitating 
isothermal equilibrium filament ( Ostriker 1964), they are consis- 
tent with some models of magnetized filaments in gravitational 
virial equilibrium ( |Fiege & Pu dritz 2000). Another attractive ex- 



planation of the observed p ~ r" z profiles is that the filaments 
are not strictly isothermal and that some of them are collapsing. 
Nakamu ra & Umemura| (p"999) have shown that the cylindrical 
version of the Larson-Penston similarity solution for the col- 
lapse of a filament has an outer density profile which approaches 
p oc r~ 2 when the equation of state is not isothermal but poly- 
tropic (P oc p?) with y^l. The dust temperature profiles derived 
for the IC5146 filaments generally show a slight, but significant 
temperature drop toward the center of the filaments (cf. Fig.|2j)), 
suggesting that a polytropic equation of state with y ^ 1 may in- 
deed be more appropriate than a simple isothermal assumption. 

The lack of anti-correlation between filament width and cen- 
tral column density (see Fig. [3]) provides another strong con- 
straint on filament evolution. If the filaments are initially formed 
with a characteristic thickness ~ 0.1 pc as a result of turbulent 
compression (see above), then the left-hand (subcritical) side of 
Fig. [3] can be readily understood. The right-hand (supercritical) 
side of Fig. [3] is more surprising since one would naively ex- 
pect an anti-correlation between width and central column den- 
sity (e.g. W oc 1 /£ ) for self-gravitating filaments contracting at 
fixed mass per unit length. However, the velocity dispersion <x v 
and the virial mass per unit length, Mii ne9V i r = 2cr 2 IG (replacing 
Mine,crit m the presence of nonthermal motions - see Fiege & 
Pudritz 2000), may increase during filament contraction as grav- 
itational energy is converted into kinetic energy. If approximate 
virial balance is maintained, the width ~ cr 2 l (G£o) ma Y thus re- 
main nearly constant as in Fig. [3] A tren d such as cr Y oc S - 5 has 
indeed been found by He yer et al. (2009 ) in their study of veloc- 
ity dispersions in Galactic molecular clouds. To account for the 
right-hand side of Fig. [3] we thus hypothesize that gravitationally 
unstable filaments may accrete additional mass from their sur- 
roundings and increase in mass per unit length while contract- 
ing. Interestingly, mass accretion onto the densest filaments of 
Taurus and Cygnus X may have been observed in the form of CO 



striations aligned with the local magnetic field ( Goldsmith et al. 
|2008l|Nakamura & Li||2008l |Schneider et al.||2010| ). Molecular 
line observations (cf. |Kramer et al.|1999| ) of a broad sample of 
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Fig. 3. Mean deconvolved width (FWHM) versus mean central 
column density for the 27 filaments of IC5146 shown in Fig [I]}, 
along with 32 filaments in Aquila and 31 filaments in Polaris, all 
analyzed in the same way (see Sect. 3). The spatial resolutions 
of the column density maps used in the analysis for the three 
regions are marked by the horizontal dotted lines. The filament 
width has a typical value of ~ 0.1 pc, regardless of central col- 
umn density. The solid line running from top left to bottom right 
shows the central (thermal) Jeans length as a function of central 
column density [Aj = c 2 / (GS )] for a gas temperature of 10 K. 



Herschel filaments in several regions would be extremely valu- 
able to confirm the validity of the scenario proposed here. 
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Appendix A: Effects of distance and dust opacity 
uncertainties and influence of the viewing angle 

There is some ambiguity concerning the distance to IC5146. The 
default distance adopted in this paper corr esponds to the valu e 



Lada et al. 



(1999) 



of 460^ pc derived from star counts by 
However, other studies placed the IC5146 cloud at a distance of 
950 ± 80 pc based on ph otometry and spec tra of late-B stars in 
the IC5146 cluster (e.g. |Harvey et al.|2008") . If the true distance 
of IC5 146 were ~ 950 pc instead of ~ 460 pc, then the widths of 
the IC5146 filaments would all be a factor of ~ 2 larger than the 
values listed in Col. 6 of Table 1, leading to a median width of z 
0.2 ± 0.06 pc instead of 0.1 ± 0.03 pc. This would also broaden 
the distribution of FWHM widths for the combined sample of 90 
filaments in IC5146, Aquila, and Polaris (see dotted histogram 
in online Fig.JTJ)). As can be seen in online Fig.[7]3, our main re- 
sult that the typical filament width is ~ 0.1 pc to within a factor 
of ~ 2 would nevertheless remain valid. We also stress that the 
distance uncertainty has no effect on the shape of the radial col- 
umn density profiles (Sect. 3.2 and Fig. [2]) or on the absence of 
an anti-correlation between filament width and central column 
density in Fig. [3] 

The following (temperature-independent) dust opacity law 
was assumed in our analysis of the Herschel data: 
k v = 0.1 x (v/1000 GHz) 2 cm 2 /g, where v denotes frequency and 
k v is the dust opacity per unit (gas + dust) mass column den- 
sity. This dust opa city law is very similar to that advocated by 
Hildebrand (1983 ) at submillimeter wavelengths, and i s consis- 
tent w i th the mean value k^o^ « 0.01 cm 2 /g derived by |Kramerj 
et al.| p003| ) in IC5146 from a detailed comparison of their 
SCUBA 850 fim and 450 /mi data with the near-infrared extinc- 
tion map of Lada et al. (1999]). However, Kramer et al. (2003 ) 
found evidence of an increase in the dust opacity /cgsoyu by a fac- 
tor of ~ 4 when the dust temperature Tj decreased from ~ 20 K 
to ~ 12 K, which they interpreted as the manifestation of dust 
grain evolution (e.g. coagulation and formation of ice mantles) 
in the cold, dense interior of the cloud. The dust temperature 
map derived from Herschel data (see Fig. [6^) suggests that Tj 
ranges from ~ 11 K to ~ 17 K in the bulk of IC5146, with the 
exception of the PDR region associated with the Cocoon Nebula 
where Tj reaches ~ 30 K. Assuming a linear increase in k v by 
a factor of 4 wh en Tg decreases fro m ~ 20 K to ~ 12 K as sug- 
gested by the | Kramer et al.| ( |2003| ) study, we estimate that the 
column density map shown in Fig. [T^ is accurate to better than 
a factor of ~ 2 over most of its extent. The possible dependence 
of k v on temperature has very little impact on our analysis of 
the filament profiles. For filament 6, for instance, Tj decreases 
by less than ~ 3 K from the exterior to the interior of the fila- 
ment (see Fig. [2])), suggesting that k v does not change by more 
than ~ 40% from large to small radii. The potential effect on the 
estimated FWHM width, W, and power-law index, p, of the fil- 
ament profiles is even smaller: W would increase, and p would 
decrease, by less than -3% and ~ 7%, respectively. 

Our Herschel observations only provide information on the 
projected column density profile £ b s (r) = 3^7 ^int(^) of any 
given filament, where i is the inclination angle to the plane of 
the sky and Si nt is the intrinsic column density profile of the 
filament. For a population of randomly oriented filaments with 
respect to the plane of the sky, the net effect is that X b s overes- 
timates Si nt by a factor < ^ > = f ~ 1.57 on average. This 
does not affect our analysis of the shape of the radial column 
density profiles (Sect. 3.2), but implies that the central column 
densities of the filaments are actually ~ 60% lower on average 
than the projected values listed in Col. 2 of Table 1. Likewise, 
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Fig. 4. Mean radial column density profile in log-log format of 
the subcritical filament labeled 14 in Fig.[lJ> (western side). Note 
that, given the relatively large dispersion of the radial profile 
along the filament (shown in yellow), the power-law behavior 
at large radii is less clear in this subcritical case than for the su- 
percritical filament shown in Fig. |2j:. 



the observed masses per unit length (Col. 12 of Table 1) tend to 
overestimate the true masses per unit length of the filaments by 
~ 60% on average. Although systematic, this inclination effect 
remains less than a factor of 2 and thus has little impact on the 
classification of observed filaments in thermally subcritical and 
thermally supercritical filaments. 
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Fig. 5. (a) Composite 3-color image of IC5146 (~ 1.6 deg 2 field) produced from our PACS/SPIRE parallel-mode data at 70, 250, 
and 500/mi. The color coding is such that red = SPIRE 500 /mi, green = SPIRE 250 /mi, blue = PACS 70 /mi. (b) SPIRE 250 /mi 
image of IC5146. 
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Fig. 6. (a) Dust temperature map derived from our SPIRE/PACS observations of IC5146. The resolution corresponds to the 
36.9" HPBW resolution of SPIRE at 500 /mi. (b) Curvelet component of the column density map of IC5146 (cf. Fig. [I]), in 
which the areas where the filaments have a mass per unit length larger than half the critical value and are thus likely gravitationally 
unstable have been highlighted in white. The positions of the 71 YSOs and 45 candidate bound prestellar cores identified with 
getsources (Men'shchik ov et aT]2010| ) in the Herschel images are overplotted as red stars and dark blue triangles, respectively. 
Candidate bound prestellar cores were selected among a larger population of starless cores (see blue triangles in Fig. [Tfr), on the 
basis of a comparison of the core masses with local values of the Jeans or Bonnor-Ebert mass (see Sect. 4.1 of Kon yves et al.|2010| 
for details). Note the excellent correspondence between the spatial distribu tion of the pres tellar cores and the supercritical filaments, 
in agreement with our earlier results in Aquila and Polaris (see Fig. 2 of Andre et al.|2010| ). 
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Table 1. Summary of derived parameters for the 27 filaments of IC5146 



Filament 


H 2 




1 d 


°T d W 


(FWHM dec 


> cr <FWHM dec ) (a) 






p w 




M line ^ 




[10 21 


cm 2 ] 




[K] 




[pc] 




[pc] 




[pc] 


[M /pc] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


1 


1.8 


0.2 


15.4 


0.6 


0.12 


0.02 


0.11 


0.09 


2.1 


0.05 


9 


2 


4.8 


4.4 


13.9 


0.7 


0.11 


0.03 


0.07 


0.09 


1.9 


0.06 


21 


3 


1.6 


0.3 


12.6 


0.5 


0.11 


0.02 


0.07 


0.05 


- 


- 


2 


4 


1.6 


0.5 


13.5 


0.5 


0.10 


0.06 


0.09 


0.07 


1.4 


0.02 


8 


5 


7.5 


4.5 


12.1 


0.6 


0.06 


0.05 


0.04 


0.04 


1.5 


0.01 


28 


6 


16 


7.6 


11.8 


0.9 


0.12 


0.04 


0.02 


0.01 


1.7 


0.04 


152 


7 


5.6 


2.4 


12.1 


0.5 


0.09 


0.04 


0.03 


0.02 


1.6 


0.03 


22 


8 


2.0 


0.7 


13.6 


0.5 


0.08 


0.03 


0.15 


0.01 


1.5 


0.05 


6 


9 


3.3 


1.3 


12.5 


0.3 


0.14 


0.02 


0.05 


0.06 


1.5 


0.04 


10 


10 


3.3 


0.6 


13.6 


0.4 


0.14 


0.03 


0.10 


0.02 


2.1 


0.06 


14 


11 


5.3 


2.5 


11.7 


0.5 


0.08 


0.03 


0.05 


0.03 


1.9 


0.04 


15 


12 


12 


4.1 


11.7 


1.0 


0.18 


0.03 


0.02 


0.01 


1.5 


0.03 


41 


13 


6.9 


5.8 


12.7 


0.8 


0.08 


0.03 


0.04 


0.09 


1.6 


0.02 


23 


14 (d) 


1.5 


0.8 


13.7 


0.5 


0.13 


0.03 


0.21 


0.01 


1.5-2.4 


0.04-0.14 


13 


15 


2.1 


2.2 


14.7 


1.1 


0.11 


0.04 


0.16 


0.06 


- 


- 


4 


16 (e) 


3.3 


0.9 


26.9 


2.1 


0.12 


0.02 


0.19 


0.05 


- 


- 


12 


17(e) 


2.0 


1.7 


18.3 


2.9 


0.09 


0.05 


0.20 


0.29 


1.3 


0.04 


6 


lg(e) 


3.9 


2.0 


22.2 


2.6 


0.11 


0.04 


0.12 


0.09 


- 


- 


49 


19 (e) 


3.6 


1.5 


26.5 


1.2 


0.09 


0.03 


0.17 


0.08 


1.9 


0.08 


1 


20 


0.57 


0.5 


15.5 


0.5 


0.08 


0.02 


0.24 


0.01 


1.5 


0.03 


3 


21 


1.5 


0.5 


14.0 


0.4 


0.13 


0.03 


0.09 


0.06 


1.7 


0.05 


6 


~ ~ (f\ 
22 {t) 


0.8 


0.1 


13.5 


0.3 


< 0.06 


0.19 


0.32 


0.05 


— 




1 


(d) 


1.4 


0.4 


13.2 


0.7 


0.14 


0.07 


0.09 


0.07 


1.5-1.8 


0.04-0.06 


7 


24 


5.6 


1.1 


12.3 


0.2 


0.07 


0.03 


0.05 


0.01 






21 


25 


2.2 


1.5 


12.7 


0.5 


0.10 


0.03 


0.09 


0.01 


1.5 


0.03 


11 


26< d > 


1.5 


0.3 


12.5 


0.4 


0.10 


0.06 


0.19 


0.04 


1.5-2.1 


0.05- 0.08 


2 


27 (f) 


2.9 


1.0 


11.5 


0.3 


<0.06 


0.03 


0.03 


0.03 


1.4 


0.02 


13 


mean 


3 


1 


13.5 


0.5 


0.10 


0.03 


0.09 


0.05 


1.6 


0.03 


11 


range 


0.6-16 


0.1-7.6 


11-27 


0.2-2.9 


0.06-0.18 


0.02-0.19 


0.02-0.32 0.01-0.29 


1.3-2.4 


0.01-0.08 


1-152 


rms 


3.4 


1.9 


4.1 


0.7 


0.03 


0.03 


0.07 


0.06 


0.3 


0.02 


30 



Notes: Col. 2, Col. 4, Col. 8: Mean projected column density, dust temperature, and Jeans length estimated at the center of each filament. 
Col. 6: Weighted mean of the distribution of deconvolved FWHM widths resulting from Gaussian fits to the radial profiles at each position 
along the filament. (a) The quoted dispersion is the standard deviation of the distribution of the measured parameter along each filament. In the 
case of ( FWHM dec > (Col. 6), this dispersion (Col. 7) is larger than the error bar plotted in Fig. [3] which corresponds to the standard deviation 
of the mean deconvolved FWHM width. (b) p and correspond to the best-fit model of the form expressed by Eq. (1) to the mean column 
density profile of each filament. The fit was weighted by l/cr(r) 2 where <x(r) is the standard deviation of the radial profile at projected radius 
r, estimated from the measured (lcr) dispersion along the filament. The same bounds (shown as vertical dashed lines in Figs. [2^,[2j:,|4]for 
filaments 6 and 14) were used to fit the profiles and to calculate the mass per unit length by integrating over the radial profile. (c ^Col. 12 
gives to the projected mass per unit length obtained by integrating over the observed column density profile of each filament. The thermally 
supercritical filaments with Mi ine > Mi ine?crit are indicated in boldface, where Mi ine)Crit = 2c 2 /G ~ 20M o /pc (for a gas temperature T ~ 12 K, 
corresponding to the median central dust temperature measured here, ignoring the temperatures of the filaments located in the PDR region). 
(d) Filaments 14, 23, and 26 have asymmetric radial profiles. A separate fit was thus performed on either side of each of these three filaments. 
The values of p and 7? flat reported in Col. 10 and Col. 1 1 correspond to the best-fit values on either side of these filaments. (e) Filaments 16, 17, 
18, and 19 are located in the PDR region called the Cocoon Nebula (see Fig.[T^ and online Fig. [5]), and their derived dust temperatures may be 
overestimated. (f) Filaments 22 and 27 have unresolved widths and are shown as upper limits corresponding to 75% of the HPBW resolution 
in Fig. [3] 
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Table 2. Median values of observed angular widths and deconvolved physical widths (a) for three samples of filaments and two 
angular resolutions 



Field 


distance [pc] 


From column density maps with 36.9' 


resolution 


From SPIRE 250 jum images with 18 


.1" resolution 






FWHM obs 


cr FWHM b s 


FWHM dec 


C" FWHMdec 


FWHM b S 


O- FWHM b s 


FWHM dec 


V FWHMdec 






["] 


["] 


[pc] 


[pc] 


["] 


n 


[pc] 


[pc] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


IC5146 


460 


59 


15 


0.10 


0.03 


44 


10 


0.09 


0.02 


Aquila 


260 


94 


21 


0.11 


0.03 


97 


31 


0.12 


0.04 


Polaris 


150 


94 


24 


0.06 


0.02 


74 


20 


0.05 


0.02 


All< b > 








0.10 


0.03 






0.09 


0.02 



Notes: Col. 3 and Col. 7: Median value of the observed FWHM angular width before deconvolution. Col. 5 and Col. 9: Median value of the 
deconvolved FWHM physical width. Col. 4, Col. 6, Col. 8, Col. 10: Dispersion (standard deviation) of the distribution of widths for each 
filament sample. (a) The DisPerSE algorithm (Sousbie 20fTJ we used to trace filaments in the column density maps derived from Herschel 
data does not consider filamentary width in its process of identifying filaments. While it is difficult to assess the completeness of our census of 
filamentary structures without dedicated tests (which will be the subject of future work), there is in principle no bias toward selecting structures 
of similar width with DisPerSE. (b) The last row refers to the combined sample of 90 filaments observed in IC5146, Aquila, and Polaris. 
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Fig. 7. (a) Distribution of deconvolved FWHM widths for the 27 filaments of IC5146 (black solid histogram, filled in orange). These 
widths have been deconvolved from the 36.9" HPBW resolution of the column density map (Fig. [T^) used to construct the radial 
profiles of the filaments. The median width is 0.10 pc and the standard deviation of the distribution is 0.03 pc. The black dotted 
histogram shows the distribution of deconvolved FWHM widths measured on the SPIRE 250 //m map which has a factor ~ 2 better 
resolution (18.1" HPBW). . This distribution is statistically indistinguishable from the distribution obtained at 36.9" resolution. The 
two down pointing arrows in the upper left mark the resolution limits for the distributions at 36.9" and 18.1", i.e., 0.08 pc and 
0.04 pc, respectively. For comparison, the blue dashed histogram represents the distribution of central Jeans lengths corresponding 
to the central column densities of the filaments, (b) Distributions of deconvolved FWHM widths for a larger sample of 90 filaments, 
combining the 27 filaments of IC5146, 32 filaments in Aquila, and 31 filaments in Polaris, all analyzed in the same way from column 
density maps with 36.9" resolution as explained in Sect. 3 and Sect. 4. The black solid histogram, filled in orange, is based on our 
default distance assumptions: 460 pc for IC5146, 260 pc for Aquila, and 150 pc for Polaris. This distribution has a median of 0.10 pc 
and a standard deviation of 0.03 pc. The horizontal double arrow at the upper left shows the range of HPBW resolution limits, going 
from ~ 0.03 pc for Polaris to ~ 0.08 pc for IC5146. The dotted histogram shows an alternate distribution of deconvolded widths 
for the same sample of filaments based on other distance assumptions: 950 pc for IC5146 (see Appendix A), 400 pc for Aquila (see 
discussion in Bontemps et al.|2010| and Appendix of Andre et al. 2010), and 150 pc for Polaris. The median value of this alternate 
distribution is 0.15 pc and the standard deviation is 0.08 pc. For comparison, the blue dashed histogram represents the distribution 
of central Jeans lengths corresponding to the central column densities of the 90 filaments (independent of distance). 



